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Radiation from a finite source in a compressible electron plasma is considered. The 
theoretical model is a perfectly conducting sphere which is excited uniformly by an annular 
slot. The configuration is such that Maxwell's equations, when combined with a (single- 
fluid) continuum theory of fluid dynamics, are separable. It is shown that a portion of 
the total power supplied by the source is radiated as an acoustic wave in the electron fluid. 

1. Introduction 

When an antenna is placed in an ionized medium such as the upper atmosphere, energy 
may be coupled into acoustic-type waves. A theory which accounts for this phenomenon 
must take account of the finite compressibility of the electron fluid. There is some evidence 
[Whale, 1963, 1964] that rocket-borne antennas excite electron-acoustic waves in the ionosphere. 
Unfortunately, the situation is complicated by the presence of the earth's magnetic field and 
the existence of a plasma sheath around the antenna. 

It is the purpose of this paper to consider an idealized model which should give some 
insight into the phenomenon of radiation in compressible plasma media. The earth's mag- 
netic field is neglected throughout. The antenna model assumed here is a perfectly conduct- 
ing sphere which is excited by a thin annular slot. The slotted-sphere antenna is assumed to 
be in intimate contact with the homogeneous plasma. 

A number of related investigations of the stated problem have been published recently. 
A short review of this work is in order. 

In a rather interesting note, Hessel and Shmoys [1962] consider the excitation of electro- 
magnetic and plasma (acoustic) waves by a point current source in a homogeneous electron 
plasma of infinite extent. The plasma is assumed to be an ideal gas coupled to the electro- 
magnetic field by the motion of the electrons. Collisions are neglected, the electron drift 
velocity is taken to be zero, and the static magnetic .field is also neglected. For these condi- 
tions, they find the ratio of total radiated power P p in the acoustic wave to that of P e in 
the electromagnetic wave to be given by 

PJP e = (1/2) (C/U) 3 (cOo/o,) 2 , 

where c is the velocity of light, u is the acoustic velocity, co is the plasma frequency, and o> 
is the operating frequency. This formula suggests that for an ionized medium, most of the 
power goes into the acoustic mode since u<CKc in most cases. This is a rather startling result. 
Hessel and Shmoys [1962] then suggest that a prescribed current distribution on a rigid sphere 
would give quite different results. 

A related problem has been treated by Chen [1963] who arrives at the same expression 
as above for the power ratio in the case of an infinitesimal electric current element immersed 
in the compressible plasma. He extends his results to a thin wire antenna of finite length. 

1 The work was supported by the Air Force Cambridge Research Laboratories, Bedford, Mass., PRO-62-201. 
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To facilitate the analysis, he assumes a sinusoidal current distribution with a propagation 
constant assumed to be identical to that for plane electromagnetic waves in the plasma. Chen 
then arrives at the remarkable result that the ratio of the powers in the electromagnetic 
and the plasma waves is independent of the velocity ratio (u/c) . 

A fundamental investigation was carried out by Seshadri [1963a] who considered a 
magnetic line source immersed in a compressible plasma. The static magnetic field was taken 
to be parallel to the infinite line source. By examining the behavior of these waves at the 
limiting case of high frequency, they are identified as the usual plasma (acoustic) and electro- 
magnetic waves. 

The excitation of coupled electromagnetic and plasma (acoustic) waves by a dipole mov- 
ing along the axis of a compressible plasma column has been considered by Yildiz and Sam- 
addar [1963]. The solution of this problem contains an acoustical shock wave which may be 
called acoustical Cerenkov radiation. 

The conversion of electromagnetic waves into acoustic waves at a vacuum-plasma inter- 
face has been considered by Hessel et al. [1962], Seshadri [1963b], Yildiz [1963], and Wait 
[1964] for various idealized models. 

2. General Equations 

The plasma medium is regarded to be a one-component electron fluid. That is, the ions 
are neglected in the equation of motion, yet their presence is required to neutralize electrically 
the plasma. It is also assumed that the amplitude of the plasma and electromagnetic oscilla- 
tions are sufficiently small that linearized equations may be valid [Oster, I960]. The average 
number density of electrons is denoted n which is regarded as constant in the plasma region. 

-> 
The pressure deviation of the electrons from the mean is p and their mean velocity is v. As 

-> -> 

usual, the electric and magnetic fields are denoted E and H, respectively. 

The linearized hydrodynamic equation of motion is 

^-> 

dv ~* 

mn ^-7=n eE—vp, (1) 

where e and m are the charge and mass of the electron, respectively. The equation of con- 
tinuity, combined with the equation of state, leads readily to 

u 2 mn V'V= — ^> (2) 

where u is the velocity of sound in the electron gas. 

Maxwell's equations for the electromagnetic fields in the plasma are given by 

vxE=_/i0 F (3) 

and 

-> 

VXH=e -—-+n o ev, (4) 

where ju and e are the magnetic permeability and dielectric constant of free space, respectively. 
Without subsequent loss of generality, all field quantities are assumed to vary as exp (i(d) and, 
thus, the derivative b/()t may be replaced by ico everywhere. 

Except where noted, collisions between particles and other dissipative effects are neglected 
in this paper. 
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3. Formulation in the Spherical Geometry 

The compressible plasma is taken to occupy the space exterior to a spherical region of 
radius a which contains the source(s). Thus, in terms of spherical coordinates (r, 6, <£), the 
tangential electric fields on the surface r=a are to be prescribed. To simplify the algebra, it is 
assumed that the prescribed field lias only an E 6 component which is independent of <j>. Thus, 
as may be verified below, the resultant magnetic field has only an H^ component everywhere. 
Furthermore, the E^ component of the field and the component of the velocity are identically 
zero everywhere. The resulting electromagnetic field may be described as transverse magnetic 
or simply TM. The continuity equation for the plasma is thus given by 

u 2 mn 



We [ siu e Vr ^)+ r i ^ sin *)]=-*^p, (5) 



while the equations of motion for the plasma are 



itornnoV r =n eE r —-^-> (6) 



io>mn Ve=n eEe— ^ (7) 



Combining (6) and (7) with (4) leads readily to the set 



Ee^-^-lirH,) e —*2 (8) 

lewr or v meoo z r od 

E r = J- -i- ^ (sin BH+) — 2 ^ (9) 

lew/* sin 00 * y meo/d/' 

— (rHJ-j— -— j^ t (10) 



m ew 2 r dr ieom mr dd 



J-^( 8in ^)-^-g, (ll) 



m eoj 2 r sin dd ieomom dr 



where 



€ 1 w . , 2 n Q e 2 

= 1 2 with cx>Z= — • 

e co me 



In the above, e/e is the relative dielectric constant of the plasma while co is the electron plasma 
frequency. 

By combining Maxwell's equations (3) and (4) with the motion and continuity equations 
(5), (6), and (7), it is not difficult to show that for r>a, 



r z? ™ +¥e [ste I (sin °z.)]+&H.=o, (12) 



and 



where 



and 



l(^)+ar^(-»if)+^=». o» 

i,=Vw»=^ (7) ' (14) 



p u 
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and 



Appropriate solutions of (12) and (13), which give rise to outgoing waves at infinity, are 

A 

H*=±^^c n P' n (co S e), (16) 

A 

P=±^0fa n P n (cose), (17) 

where c n and a n are unknown coefficients, and where the summation is over integral values 

A 
of n. In the above, K n (z) is a spherical Bessel function which satisfies 

4P^H"^]tM, OS) 

while P n (cos 6) is a Legendre polynomial which satisfies 

sin Jg (sin 6 ^)+w(n+l) sin 2 0P„ = O, (19) 



and, finally, 



PUcosd)=~P n (cosd). (20) 



Following Schelkunoff [1943, 1952], the spherical function which satisfies (18) is taken to be 

t> f n 2 ^ (n+m)! 

Thus, for k e r and ^/'>>1, 

K n {%k e r)»~e- ik *\ (22) 

and 

K n (ik p r)^e- ik p r , (23) 

which illustrates the outgoing character of the waves. It is also evident from the definitions 
of k e and k p , given by (14) and (15), that radiation will not take place if €/e <^0 (i.e., co< co ). 
In what follows, it is assumed that k e and k p are both real, corresponding to e/e >0. 

An explicit general expression for the r component of the velocity v is obtained by using 
(11), (16), and (17), and getting an assist from (19). Thus, for r^a 

v-f^fc ™(n+l) K n (ik e r) e k p V K' n (ik v r) K n (ik p r) -]\ 

Vr ~hi\ n eo> 2 mr (ik e r) a » eam rn [__ ik 9 r (ik p r) 2 jfXWOsV), (24) 

where the prime indicates a derivative with respect to the argument of the Bessel function 

A A 

[i.e., K' n (z)=dK n (z)ldz\. The boundary condition imposed on the velocity is that v r =0 at 
r=a. In other words, the spherical boundary is assumed to be rigid as far as the acoustic 
wave is considered. A simple calculation then leads to the result that 

A 

r K n (ik e a) l 

g »— -to in en * L (Hc e a) J / n 

7- n{n+1) ^k P Z*— " a ,., _/ (25) 



' V K n ' (ik p a) K n jik p a) ~j 
L (ik p a) (ik P a) 2 J 



The problem is formally solved when a n is determined by the specified conditions at the 
source. This task is deferred until later. 
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4. Radiated Power Considerations 

The total power P, radiated from the source, may be obtained by integrating the radial 
component of the Poynting vector over an infinitely large sphere surrounding the source. 
Thus, 



P= Lim i f * [EeH%+pv*]2wr 2 sin 

r->oo Z Jo 



ede. 



(26) 



The asymptotic form of the field components, valid for & e r>>l, are obtained from (8), (9), 
(16), (17), and (22). Thus, 



^ZX~P;(cos0), 
n =i tti e f 



and 



-ik p r 



Ee^-TJe X! On -JT— P^(C0S 0) , 
n = l btL e i 



(27) 



(28) 



where t\ e =\x^\k e . In a similar manner, by making use of (10), (11), (16), (17), and (23), it 
is found that 



P~ZX^r^»(cos0), 



n=l %/C v V 



and 



V r ^- 



e i 



i/« ' a 



^2a n - 7 r—P n (cos 6). 



€oM m^i n ik p r 



(29) 



(30) 



The 6 integration in (26) may now be readily carried out if use is made of the orthogonality 
relations [Morse and Feshbach, 1953] 



P n (cos 0)P m (cos 0) sin 



Odd 



if n= m 



and 



'p;(cos0)P;(cos0) si 



2n+l 
=0 iin^m, 

sin vav= V ' y if n=m 
2n-fl 

=0 ifn^m. 



Thus, it is found readily that 



P_9^ !k ^ n ( 7l + 1 ) u 12,0 *o ^ Kl 2 



(31) 



(32) 



(33) 



where the first summation corresponds to the power in the electromagnetic wave and the second 
is the power in the acoustic wave. In order to display some important properties, (33) is 
rewritten in the form 



P _r> Vb f-> n(n+l) , , 



where 



H££(2n+1) 



n(n-\-\) n m k p 



(34) 
(35) 



Using (25) and noting that o)l=n e 2 l(me ), it follows that 



A n =n(n+1) ^ 



\K n (ik e a)\ 



co 2 (& e a) (i p a) 



K n (ik v a)- 



K n (ik p a) 



which is a rather fundamental result. 



(ik p a) 



(36) 
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It is evident from the above development that A n is the ratio of the acoustic power in mode 
n to the electromagnetic power in mode n. The result is applicable to any radiating source 
which produces azimuth ally symmetric TM modes provided that the boundary at r=a is rigid 
in the sense described. It is rather interesting to observe that the ratio A„ depends on three 
parameters, namely, (i) k e a, which is the circumference of the rigid sphere in electromagnetic 
wavelengths, (ii) k p a, which is the circumference of the rigid sphere in acoustic wavelengths, 
and (iii) a>o/co 2 , the normalized plasma frequency, which is proportional to the electron density. 
It is worth noting that the first two parameters may be written 



(i 



and 



i) k 6 a=(ka)(l-^j, (37) 



(ii) k P a=(ka) (l-^J (£j> (38) 



where ka is the circumference of the rigid sphere in free-space wavelengths, u is the acoustic 
velocity, and c is the velocity of electromagnetic waves in free space. 

It is important to note that, everything else considered equal, A n is proportional to n(n-\-l). 
Thus, the higher order modes excite appreciably more acoustic-type radiation than the funda- 
mental mode (i.e., n=l). 

Some limiting cases are worth discussing briefly. If both fc e a>>l and k p a^>^>\, it 
readily follows from (22) and (23) that 

A^^+D^g)^. (39) 

Since in most practical situations u/c is of the order of 10" 3 or 10~ 4 , this quantity is very small, 
indeed. 

When k p a^>^>l, without any further restrictions, it is seen that 

which is also small unless (k e a) itself is somewhat less than unity. For example, if k e a<^<^l, 
while keeping &„a>>l, it follows from (21) that 

, , ,, a%u i n»n 2 , N 

which may not be small, particularly for the higher order modes. 

Another limiting case which is mostly of academic interest is when both (& e a)<<l and 
(k p a)«l). Then 

*■"*+! J W ' (42) 

Since ?//c<C<0, this quantity is enormous even when co/co is somewhat less than unity. Setting 
n=l for the fundamental mode, the resulting formula for A! is identical to that derived by 
Hessel and Shmoys [1962] for a point current element in a compressible plasma. From this 
particular limiting case, one may conclude that the formula of Hessel and Shmoys is valid 
only for an actual current element if its length is small compared with the wavelength of the 
acoustic-type wave. 

5. Excitation Considerations 

In the foregoing discussion the precise manner in which the rigid sphere is to be excited 
was not specified. Actually, the results derived and the conclusions drawn require only that 
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the spherical surface bounding the source(s) is effectively rigid. It is also assumed, of course, 
that only azimuthally symmetric TM modes are to be excited. One possibility would be a rigid 
dielectric sphere with a radially oriented electric dipole contained within [Wait, 1964]. An- 
other suitable method of excitation, which is considered in (his paper, is to assume that the 
source is a perfectly conducting rigid sphere of radius a. At the colatitude angle O , a uniform 
voltage Vq is impressed across an annular slot of infinitesimal width. In other words, the 
tangential electric field on the spherical surface is specified by the condition 



E t (e,a)=(y l ja)s(e-e ), 



(43) 



for r=a and O<0<7r, where 5(0 — O ) is the unit impulse function at 6=6 . The correctness 
of this result is confirmed by noting that 



f 



Ee(d, a) add =Lim 



f 0o+ a 

im 

->0 J do- A 



E B (6, a)add=Vo } 



(44) 



where 2A is the width of the slot which is allowed to approach zero. 

A general expression for the E% field in the exterior region is obtained by using (8), (16), 
and (17). Thus. 



A 



e K n (ik p r) 



n^b L (ik e r) eu 2 mr (ik p r) 

In the limit as r— »a, this result may be expressed in the form 



0»] P'n 



(cos 6). 



(45) 



where 



with 



tie(a,d) = ^2A n P' n (cosd) 1 

n = l 



A 

K' n {ik e a) 

A "= —ye „-!. ~ C„{1— 5 n ), 



k v ic7na 



ik e a 

A 

e /a n \ K„(ik p a) 






K'„ (ik e a) 
Using (25), an explicit formula 8„ is obtained which may be written 



(46) 



(47) 



(48) 



<* (^a)(M t n (ik e a) 



A 

K n (ik p a) 



tdha)- 1 ^^ 



ik v a 



(49) 



which is valid for the rigid boundary condition at r=a. 

By making use of the orthogonality property of the Legendre polynomials, as indicated 
by (32), the coefficient A n , from (46), may be written 



A »=2^+l) Jo Eeia ' 6)K(C0S 6) Si " M8 - 
Then, as a consequence of (43), 

An 



271+1 V D , , Q \ ■ a 

"2^+1) T Pn(c0S ^ )sm *" 
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(50) 



(51) 



Explicit expressions for the electromagnetic field components may now be obtained if (47) 
and (51) are combined with (16), (17), and (25). For example, 

jr _ Vp sin fl » (2n+l) K n (ik e r) P' n (cos 6> ) 

H - 5T gig+ij ^ -r=ar ^^^ <52) 

is an exact formula for the magnetic field which has only a <t> component. In the far zone, 
where \k e r\^>^>l, this simplifies slightly to 

H _ VpsmOoe-**'' « (2n+l) P;(cos d )P' n (cos 6) ^ 

(l—8 n )K n {ik e a) 

which is a spherical wave behaving like e~ ikeT/r . In the far zone, the 6 component of the electric 
field is obtained from 

Ee^rjeHj where 7 le =» a)/k e =120<jr (l ° ) • (54) 

The total power P e in the electromagnetic wave is obtained from 

P e =\ Re Lim P ' E e H*(2>Kr 2 ) sin $d9, (55) 

" r->°o Jo 

which is evaluated to give 

Pe=7i — sim v 2_j — ( — tt\ 1 * oo) 

2,. £i»(n+l) i!^,,^^),, 

Referring back to (34), it is evident that the total power radiated, P (electromagnetic and 
acoustic) , may be written as the sum 

P=P e +P P , (57) 

where the radiated power in the acoustic wave is given by 

l Ve -iM»+lJ \l-8 n \'\(K' n (iha)\* 

where A n is given by (36). 

As an important check on this derivation for the total power P, the Poynting vector may also 
be integrated over the area of the annular slot. This latter method which is related to the 
so-called "emf method" in antenna theory [e.g., Wait, 1959], is not very rigorous but it does lead 
to correct results if applied with care. Therefore, the power P should also be obtained directly 
from 



P 2 

Area of slot 

This may be written 



\ f TLe[E e H*] r=a da. (59) 



P= I Lim P° +A Re [E e H*] r=a (2Ta 2 sin 6 )dd (60) 

=wa sin e Q V Re (fl*) r - fl , (61) 

where V is the voltage across the slot (which may be regarded as real). Using (52), it is seen 
that according to this approach 

A 

D _ 7T sin 2 O 



|^ VI ± ^ [P'Acos ,,]■ Re |" A K ^ a) ( 62 ) 

2,. „.m(n+l) lk'MMlSn)_ 
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The exact correspondence between (57) and (62) is a result of the identity 



Ke 



~K n {ik e a) 1 



A 1 —8 

K' n (ik e a) " 



1+A„ 



\K' n {ik e a)\ 2 \l-t> n 



(63) 



where A n and 8 n are defined by (36) and (49), respectively. This identity is verified without 
difficulty, if it is remembered that, 

Rek n (iz)[K' n (iz)]*=-l, (64) 

for z real. 

Equation (62) may be used for the power calculation for a slotted sphere of arbitrary size 
immersed in a compressible plasma. When k e a and k p a are not large compared with unity, 
the series converges rapidly. Then, it is of interest to see that the power in the first mode may 
be written 



where, according to (36), 



and according to (49), 



3 7rsin 4 fl 1+Aj (J y 



A x ^2 



C0 



(Kay 



(k v a) 



K[(ik v <i) K ^ a) 



(ik p a) 



* . <4 2 

w 2 (k p a) 



A 

Ki(ik v a) 



km P a) ksk ^ 



ik p a 



(65) 
(66) 

(67) 



where it should be remembered that 



k v a=(c/u)k e a. 



(68) 



If u/c is sufficiently large that k v cb^>^>l, while still maintaining & 6 a<^<0, the above formulas 
for Ai and 5i simplify to 



W^-- 



and 



co 2 c (k e aY 



«,s-2i=| 



COq U 1 



c k e a 



(69) 



(70) 



Using the above result, it is possible to write a simple formula for the ratio of the total power P 
for a compressible plasma when u/c is finite to the total power P 9 for an incompressible plasma 
(such that uje=0). Thus, by using (65), (69), and (70), it is seen that, for n=l, 



P 



1+2^- 



c {k e a,y 



<4 fu\ 2 _J__ 
<o 4 W (k e ay 



(71) 



1+4 



For this same condition, it is worth noting that the ratio of the power P e of the electromagnetic 
wave in the compressible plasma to the total power for the incompressible plasma, is given by 



P.. 

Po = 



o) 4 \c/ (k e a 
1135 



(72) 



As indicated previously, the ratio of the acoustic power P p to the power electromagnetic 
P e , for n=l, is 

^2 4-^ (73) 

P e u 2 c (k e ay 

for these same conditions (i.e., k e a<^<^\ and k p a^>^>l). 

6. Concluding Remarks 

On the basis of the analysis in this paper, it would appear that an electron-acoustic wave 
will be excited by a slotted-sphere antenna which is immersed in a compressible plasma. 
However, if this effect is to be measurable, the operating wave frequency must be an 
appreciable fraction of the plasma frequency. It also appears that the relative power in 
the acoustic wave is increased as the overall dimensions of the antenna are reduced. 
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